Experimental realization of a fetching algorithm in a 7 qubit NMR quantum computer 
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Searching for marked items from an unsorted database is an important scientific problem and a 
benchmark for computing devices as well. Using a 7-qubit liquid NMR quantum computer, we have 
demonstrated successfully an hybrid quantum fetching algorithm that finds marked items using only 
a single query. The essential idea is the operation of quantum computers in parallel. We gave the 
detailed pulse sequence for coherent control of the 7 qubits. The pulse sequence demonstrated here 
is not only useful for ensemble quantum computation, but also can be regarded as a general purpose 
control-gate which is useful for experimental design of quantum algorithms and general quantum 
information processing task in other quantum computer schemes. A generalization of the algorithm 
that is scalable to arbitrary qubit number is also provided. 



The perspective enormous power of quantum computers(QCs) has sparked intensive efforts worldwide. Various 
schemes have been put forward. Liquid NMR quantum dot cavity quantum electrodynamics(QED) 

Hnear ion trap 0, superconducting qubit solid-state NMR are just a few examples in this long list. By the 

way they measure, they can be put into two main categories: single qubit-measurement quantum computer (SQC), 
and ensemble qubit-measurement quantum computer(EQC). Typical examples are Kane's proposed silicon quantum 
computer(SQC), and the liquid NMR(EQC). SQC is the ideal system for quantum computation. EQC is in fact 
an ensemble of many SQC copies. The control and measurement in EQC, such as liquid NMR, is easier than SQC 
for existing technology. Thus it advances farther in demonstrating quantum information processing such as Dcutsch- 
Josza [p^-p^, quantum search and counting 1 6 - 19| , o rder-finding pof, q uantum Fourier transformatio n pT 22 1, Shor's 
factorization , quantum error correction code |24|] , dense coding [p5[ , and GHZ-state preparation [p6[. 

Usually an ensemble is in a mixed state. Effective pure state(EPS) technique is employed to prepare the 

system so that the ensemble behaves like a pure state. In liquid NMR, a measure that indicates to what extent the 
system is close to a pure state is the polarization e [Q. If the polarization is equal to 1, the ensemble becomes a pure 
state ensemble, EQC is equivalent to SQC. By increasing the magnetic field or employing polarization transfer p7,^, 
polarization can be increased. At present, NMR-QC is operating at a regime with e ~ 10^^. Liquid NMR QC may 
be used to investigate the behaviour on representative systems and study problems that will be encountered in future 
scalable technologies. 

Redundancy in EQC not only makes qubit control and measurement easy , but also provides source of additional 
computing power. Madi, Briischweiler and Ernst p9[ | put forward Liouville space quantum computation(L-QC). In 
this mode, the quantum parallelism and the classical parallelism in an ensemble are exploited, and it can achieve 
exponential speedup for certain problems which are impossible in a pure quantum computer. Different from EPS 
where all the sub-ensembles are made into the same state as much as possible, L-QC exploits the diversity of states of 
constituent molecules in an ensemble. Notably, unsorted database search can be exponentially fast in L-QC. Unsorted 
database search is a benchmark problem for computing devices. It is a mathematical NP-hard problem and plays 
an important role in computational complexity theory. It is widely used in science and technology, for instance it 
can break DES-like encryption schemes ]3l| ]. To find a marked item from a database with A^-items, it requires about 
A^/2 queries in a classical computer. In a quantum computer, Grover's algorithm [ p2[ uses 0{\/N) steps, though a 
significant square-root speedup but still exponentially slow. When searching a database with unknown number of 
marked items, an estimate of the number of marked states is required because it depends sensitively on this number. 
Furthermore it has been shown that Grover's algorithm is optimal for quantum computers ^0 33 3^ . This optimality 
restriction can only be broken off from outside of quantum computer and it has been achieved by Briischweiler in 
L-QC Briischweiler 's algorithm requires only log2 A^ search steps. For instance, to decipher a DES code, it is 

equivalent to find a marked item from a database with A^ = 2^^ ~ 10^^ items. A quantum computer running Grover's 
algorithm requires about 0.8 x 2^^ ~ 2 x 10^ steps, whereas Briischweiler's algorithm requires only 56 search steps. 
Recently, it is found that L-QC computing can achieve the absolute maximum in unsorted database search, a single 
query suffices to find every marked item using a hybrid quantum fetching algorithm | ]36[ | . 

In this paper, we report the experimental implementation of this hybrid quantum fetching algorithm in a NMR 
system with 7 qubits. The implementation requires coherent quantum control of all the 7 qubits. Pulse sequences for 
the quantum gate operations are given explicitly. We also scale up the algorithm to arbitrary qubit number which 
can be used in scalable ensemble quantum computer. In particular, we found that the algorithm can be implemented 
directly from thermal equilibrium. 

Here we briefly introduce the fetching algorithm, details can be found in Ref. |36[. In operator-product formalism 
, a pure-state \4>m) is represented by direct product of the polarization operators |29|, \4>m) = \aap...af3) ■^=^ pm = 
/f where |a), are the spin up and down eigenstates of Iz, the Pauli spin operator respectively. 



The one-qubit state |a) and \(3) can be expressed in terms of the polarization operator, namely 2/" = 2|Q!)(a| = 
(1 + 24), 2lP = 2|/3)(/3| = (1 - 2h). In L-QC with 7-qubit, one qubit is used as ancilla as required by L-QC, the 
other six qubits represent a database with 2^ = 64 items, namely |0) = |000000), |1) = |000001), ...,|63) = |111111). 
To start, the L-QC system is prepared in a mixed state, apart from a scaling factor 

po = /o" = |0, 0) (0, 0| + |0, 1) (0, 1| + ... + |0, i) (0, z| + ... + |0, 63) (0, 63|, (1) 

where we have written the ancilla qubit state explicitly for clarity, namely jOOOOOOO) — |0, 0). Then the query function 
/ is applied to the system, which changes state into 

Pi = 1/(0), 0) (/(O), 0| + ... + |/(z), z) {f{€),i\ + ... + 1/(63), 63) (/(63), 63|. (2) 

For constituent (to distinguish pure state from mixed state, we call a part in a pure state as component, and a part in 
a mixed state as constituent) corresponding basis i, query function does nothing to it if i does not satisfy the query, 
otherwise it flips the ancilla qubit. After the implementation of query, a read-out pulse is applied to the ancilla qubit 
and its spectrum is measured. Those constituents that satisfy the query have their ancilla qubit flipped, and on the 
ancilla qubit spectrum, the peaks corresponding to these constituents will have their peaks inverted. The item that a 
peak corresponds is determined by the peak frequency. Then by looking at the ancilla qubit spectrum, one can easily 
read out the marked items. 

The experiment is carried out on a VARIANT INOVA-600MHz spectrometer. The liquid sample is ^'^C labeled 
crotonic acid with formula C^H^C^H^ = H^C'^02H , which is manufactured by Cambridge Isotope Laboratories 
Inc. The same substance has been used in Ref. Q to prepare GHZ state. The chemical structure and the J couplings 
can be found in Fig. 2 of Ref. |Q. The seven qubits are the four carbon nuclear spins, two proton spins in the middle 
connected with the middle two carbon. The 3 methyl proton spins(denoted M) have the same chemical shift, and 
they are treated as a single qubit. The sample is firstly dissolved in a 0.5 ml acetone and then injected into a 5 
mm tube, deaired and sealed. The ancilla qubit is chosen as C^. Nuclear spins of H^, C^, C^, H^, and are 
assigned as qubit 1, 2, 3, 4, 5 and 6 respectively, in decreasing order by the J magnitude of respective nuclear spin 
with the ancilla bit. Our measured J values agree very well with that in Ref. p^ . The Jgi coupling between qubit 
and the rest six qubits are(in Hz): 156.0, 69.7, 41.6, —7.1, 1.4, —0.7 for i from 1 to 6 respectively. If J^i is positive 
such as qubit 1, 2, 3 and 5, we assign \a) = |0) and \(3) = |1), and if Joi is negative such as qubit 4 and 6, we assign 
|a) = |1) and = |0). In doing so, the peaks in the ancilla qubit spectrum will represent numbers monotonically 
in increasing order starting from jOOOOOO) with frequency ujq -f- "'I'^ojl on the far left to 1 111111) with frequency 

^0 - Z^Li ■^I'^oil on the far right. 

The query function / is important in unsortcd database search. It is treated as a blackbox ||,||J3|. The query 
blackbox used in this work is given in Fig.|l|a. It checks whether an items's first 3 bits values are 100. The H-gate in 
the figure is a Hadmard-like gate with the expression 

TT « / 1 1 



H° = exp(-z IT II) exp(-z - J^;) = _ } j . (3) 



It is realized using NMR pulse sequence B!^{'k)R^{^). The superscript (subscript) is the qubit (direction) on which the 
pulse is applied. The triple-qubit controlled phase rotation in the query is explicitly as 

«rM^»p(-„/«<i^(i^ii^,. ,4) 

M_2/i) Q^2/^) 

It is clear that when the first three qubits are in values 1, and respectively, the expression -^^ — y^'^ — T^'^ — 
gives an identity operator, and the whole expression in (Q) becomes a rotation about the z-axis through tt. Otherwise, 
one or more of the factors becomes zero, and Eq.(|^) reduces to an identity operator. This implements the triple-qubit 
controlled phase rotation. Further it can be reduced to basic one- and two-qubit gates by the recursive formula 

exp(-iA2"/fe,,4,,.../fe„,/fe„^,,) = K exp(-zA2"-i/fc,,.../fe„_,,/fe„^,,)l/+, {n > 2) (5) 

given in Refs. ||3§], where 



TT TT TT 

Vn = cxp(-i-/fe„^j^) exp{-iTTlk^^Ik„_^^^_) exp(i-/fc„^j:r) exp(-j-/fe„^jj^). (6) 



The whole pulse sequence is drawn in Fig.|| It involves the coherent control of all the 7 qubits. It consists of 47 
spin-selective radio frequency pulses, together with several free evolution intervals. The decoupling is implemented 
using the compound decoupling pulse sequence WALTZ 16 Eo|. The total duration of the pulse sequence is about 90 
ms which is well within the system's 2s decoherence time (7^ [ p6| . 

We actually implemented the algorithm directly on the thermal equilibrium state without preparing 1^ . In thermal 
equilibrium, the density operator is 



(7) 



The state we need is Iq = {1 + cr^) ~ (7° and can be written as p' = (I — /37o-Bcr^), apart from an identity operator 
and a scaling factor. 7^ is the gyromagnetic ratio of qubit i, and /3 — 1/kT. After the implementation of query, only 
the ancilla quibit state is changed, 

Ufp'U] = ^ - PUn^alu} - /3 1] i.Bal^ . (8) 

Upon measurement, a selective read-out pulse is applied to the ancilla qubit, the first and last terms in eq.(H) do not 
contribute to the signal. The signal we acquire is only from the second term which is the one we need. Thus we can 
perform the algorithm directly from thermal equilibrium. Shown at the top of Fig.^ is the ancilla qubit spectrum at 
the thermal equilibrium state. It represents the unsorted database. The spectrum is divided into 8 major groups, 
and each group contains 4 peaks where two inner peaks are about 3 times as high as the outer two short peaks. 
This is because qubit 4 is composed of 3 identical proton spins, and they produce 4 energy levels and thus make 4 
modulations on the ancilla bit spectrum with modulation frequency 3Jo4 (all three spins up), J04 (two spins up and 
one spin down, triply degenerate), — Jo4(two spins down and one spin up, triply degenerate), and — 4Jo4 (all three 
spins down). Because the multiplicity of levels due to qubit 4, the number of peaks in spectrum is doubled. Thus 
within each group the left two peaks(one short and one tall) corresponds to qubit 4 at state 0, and the right two peaks 
corresponds to qubit 4 at state 1. While the doubling of levels gives additional work in preparing an effective pure 
state in EPS QC, it is advantageous for L-QC because this peak-doubling makes the identification of peak inversion 
a lot easier. Within each peak in Fig.^, there are still 4 finer peaks due to the different states of qubit 5 and 6. In 
finer frequency resolution, they are clearly separated and visible. All the 128 peaks in the spectrum correspond items 
from at the far left and 63 at the far right. It is interesting to notice the separation of the major groups. Counting 
from left, the distance between group 2 and 3, 6 and 7 are shorter than that between groups 1 and 2, 3 and 4, 5 and 
6, 7 and 8. The distance between 4 and 5 is the largest. These distance differences can be easily understood from the 
J parameter magnitude. The centroid positions of these peaks are of the form ±|Joi| i | J02I i | Jos]. 

Shown in bottom of Fig.H is the ancilla bit spectrum after the implementation of the query. The spectrum is a 
result with only 8 scans. Shown in the middle of Figjl is part of the spectrum after the query, here it is clearly shown 
that there are 4 peaks within each peak. Even with this small scan number, it is clearly visible from the figure that 
the peaks in the fifth group have their peaks inverted. The doubling of peaks due to qubit 4 makes the identification 
of peak inversion easier, in particular by looking at the inner tall peaks. They correspond to items from 100000 
to 100111 which are the items that we are searching. As compared with that before the query, the height of peak 
is reduced. This is due to the decoherence that occurs during the computation. The result is striking because no 
measures have taken to compensate the decoherence. The fetching of the marked states is successfully demonstrated. 

It is worth pointing that the algorithm demonstrated here is in fact the parallel operation of quantum computers. 
The capability of doing quantum computing in parallel lies in the unitarity nature of QC operations. A function / in 
a QC is carried out by a unitary transformation Uf which corresponds to a sequence of radio frequency pulses and free 
evolutions. This pulse sequence evaluates the function value for |00...0), and the same pulse sequence also evaluates 
the function value of |ll...l) and any other inputs! When applied to a superposition of basis states, it evaluates the 
function values of all the basis states! When applied to an ensemble in mixed state, the pulse sequence transforms each 
constituent sub-ensemble of basis states into their corresponding function values. The pulse sequence is ubiquitous to 
all molecules in the ensemble and acts as a super-commander to all the basis states. Each basis state goes through 
the transformation required by the specific function irrespective of its surroundings: alone, or lined up hand by hand 
with other basis states in a superposed pure state, or scattered here and there in different sub-ensembles in a mixed 
state. 

The algorithm can be scaled up to systems with arbitrary n qubit system. An example for the query network for 
n-qubit system is shown in Fig.Qa. The multiple-qubit(iiZ2..in) controlled phase rotation can be written explicitly 

i?f (^) = exp(-^|^/°(l + s,{~\y^2ll)...{l -f s„_i(-l)^"-2/r^)(l + s„(-l)'"24")), (9) 

where Si denotes the sign of Joi. Using formulas (^) and (H), it can be decomposed into basic one- and two- qubit 
gate operations. They might be used in other ensemble quantum computers, such as the proposed schemes in Refs. 

00- 

To summarize, we have implemented a hybrid quantum fetching algorithm in a 7 qubit NMR quantum computer. 
The algorithm is essentially the parallel operation of quantum computers. A scale-up version of the algorithm is 
provided. Furthermore it is demonstrated that the algorithm can be implemented in a thermal equilibrium state 
directly. 
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FIG. 1. Quantum network for query: (a)used in this work; (b) in arbitrary qubit number system. _ff is a Hadmard-like 
transformation. A solid dot at the intersection with a qubit line means 1, and circle means 0. Rz(pi) is a tt rotation about the 
z-ajds. 












4 


T 








1 




1 
1 




.1 
1 








2 




- 1/2J,„ — 




— 1/2J» - 
1 1 

1 1 


— mho - 




— 1/2J50 H 
1 1 

1 1 


— lOJio — 




— 1/2J3() — 




— 1/2J21, — 




— 1/2J5,, 




1 


3 










T-L-L 






















1 








— ^ 






— L- 










— L 








5 








waltz 16 


waltz 16 




waltz 16 






waltz 16 




waltz 16 




waltz 16 




6 








waltz 16 


waltz 1 6 




waltz 16 






waltz 16 




waltz 1 6 




waltz 16 




4 








waltz 16 


waltz 16 




waltz 16 






waltz 16 




waltz 16 




waltz 16 





FIG. 2. The pulse sequence for the fetching algorithm. The symbols(in brackets with sub-total number) for different pulses 
(in degrees) are: 180 along j/(thick long line, 23), 90 along y(thin short line with a bar at top,l), 90 along — j/(thin short line 
with a bar at middle, 7), 112.5 along j/(thin long line with bar at top, 4), 22.5 along — i/(thin long line with bar at middle, 1), 
67.5 along y{thm long line with arrow, 3) and 90 along a;(thin short line, 8). 
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FIG. 3. The spectra before(top) and after(bottom) the implementation of query. The plot in the middle is a part of the 
spectrum boxed at the bottom. Here it is seen that each peak in the bottom spectrum has 4 sub-peaks. 



